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We show that the order-disorder phase transition in the three state Potts ferromagnet on a square 
lattice is driven by a coupled proliferation of domain walls and vortices. Raising the vortex core 
energy above a threshold value decouples the proliferation and splits the transition into two. The 
phase between the two transitions exhibits an emergent U(l) symmetry and quasi long range order. 
Lowering the core energy below a threshold value also splits the order-disorder transition but the 
system forms a vortex lattice in the intermediate phase. 


Phase transitions in a variety of systems are driven by 
the proliferation of topological defects m- Manipula¬ 
tion of defects, therefore, provides a natural route to¬ 
wards altering the nature and location of phase tran¬ 
sitions, which in turn can significantly alter the phase 
diagram itself. The role played by a single type of de¬ 
fect, and the effect of manipulating it, has been studied 
in models of superfluids Ml- liquid crystals HI! and 
Heisenberg ferromagnets [1QHT2] . In a large class of sys¬ 
tems, however, the phase diagram is determined by the 
proliferation of not one but multiple types of defects 
[20] . We would like to identify a minimal model in which 
the interplay between two types of defects, and the effect 
of manipulating them, can be studied clearly. 

The two state (Ising) ferromagnet on a square lattice 
is one of the simplest spin models which exhibit a de¬ 
fect driven phase transition. Domain wall defects ap¬ 
pear as small loops in the ordered phase of the model 
and drive a transition to the disordered phase upon pro¬ 
liferation HU [22] . The next simplest model, the three 
state Potts ferromagnet, also exhibits an order-disorder 
transition [23[ [24]. This model, however, supports the 
formation of domain wall as well as Z 3 vortex defects. 
An approximate energy versus entropy balance calcula¬ 
tion suggests that the system disorders because vortex- 
antivortex pairs unbind as soon as the domain walls pro¬ 
liferate HU. Apart from this calculation, the role played 
by the two types of defects in the model has remained 
largely unexplored. 

In this Letter, we use Monte Carlo simulations to 
demonstrate that the interplay between domain walls and 
vortices in the three state Potts model generates a rich 
phase diagram [Fig. [I]. We show that the order-disorder 
transition in the model is driven by a coupled prolifer¬ 
ation of the two types of defects. When we raise the 
core energy of the vortices by an amount A, the model 
continues to exhibit the order-disorder transition upto a 
certain threshold A = A + . Above A+, the vortices prolif¬ 
erate after the domain walls and the transition splits into 
two. The phase, which appears intermediate between the 
two transitions, exhibits an emergent 1/(1) symmetry and 
quasi long range order (QLRO). When we lower the core 
energy using negative values of A, the order-disorder tran¬ 
sition of the pure Potts case becomes sharper. Below a 



FIG. 1. (Color online) Schematic phase diagram in the 
parameter space of temperature T and vortex suppression A. 
The pure Potts model corresponds to A = 0. Estimates for 
A+ and A_ are given in the text. Domain walls (black lines), 
vortices (blue dots) and antivortices (red dots) are overlaid 
on typical spin configurations obtained in each phase. 


threshold A_, the transition again splits into two. In this 
case, however, the intermediate phase is a vortex lattice 
in which the vortices and antivortices display sublattice 
ordering. 

Before discussing the phase diagram in detail, we de¬ 
scribe how the defects are identified for a given config¬ 
uration of spins. Each spin <7^, at vertex i on a square 
lattice A, can be in one of three states: <7^ G {0,1, 2}. Do¬ 
main walls and vortices reside on the edges and vertices, 
respectively, of the dual lattice A' which is a square lat¬ 
tice displaced from A by half a lattice spacing along each 
axis. If two spins across an edge (i, j) G A are in dissimi¬ 
lar states, then a domain wall is placed on the dual edge 
in A'. The vorticity at each dual vertex i' G A' is deter¬ 
mined by calculating a discrete winding number ujp na 
For Z 3 vortices, cup = (A ba + A cb + A dc + A ad )/3 where 
A ba represents (a b — a a ) wrapped to lie in [—1, +1] and 
<7 a , a b , <7 C , (j d are the four spins on the square plaquette 
in A surrounding i’ in an anticlockwise sense. A vortex 
(antivortex) is present at i' if uv is +1 (—1). 

A standard method for suppressing the formation of 
vortices in models of superfluids involves raising the vor¬ 
tex core energy by an amount A 0 0 H3 ng. Upon 
inclusion of such a term for Z3 vortices, the three state 
Potts Hamiltonian, with nearest neighbor ferromagnetic 
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FIG. 2. (Color online) The order-disorder transition of the pure Potts model (A = 0), marked by the decay of S and divergence 
of xs, is accompanied by a simultaneous increase of pd w and p vx . Weak suppression of vortices (A = 2) shifts the transition 
to a higher temperature. Stronger suppression decouples the simultaneous proliferation of the defects and splits the transition 
into two. The data has been obtained for L — 16 (diamonds), L — 32 (triangles) and L — 64 (circles). 


interaction J > 0, becomes 

n = j (i-<W + a E M- (i) 

(i,j)E A i'G A' 

If the number of domain walls and vortex defects cor¬ 
responding to a given spin configuration is denoted by 
Ndw and N VX: respectively, then 0 can be rewritten as 
1~L = JNd w + \N VX , clearly indicating that the statisti¬ 
cal behavior of the model depends solely on the number 
of defects. In particular, the behavior depends on the 
number of domain walls and their three branch intersec¬ 
tions. Each i' G A' is visited by zero, two, three or four 
domain walls. Out of these four scenarios, ujp ^ 0 only 
when three domain walls intersect. The present model 
can, therefore, be equivalently expressed as a domain wall 
loop model m with a fugacity parameter A controlling 
the density of three branch intersections. 

We have determined the phase diagram of the model in 
the two-dimensional parameter space of A and tempera¬ 
ture T by simulating 0 on a Lx L square lattice. As the 
plaquette based A term cannot be incorporated into cur¬ 
rently known cluster algorithms, the spins were updated 
using a single spin-flip algorithm [28] . We measured the 
density of domain walls p^w = ^dw/^L 2 , the density of 
vortices p vx = N vx /L 2 and the Potts order parameter 
S = 3(max{no, ni, 77 - 2 } — 1/3)/2, where n a represents the 
fraction of spins in state a. Large autocorrelation times, 
arising from the use of the single spin-flip algorithm, were 
estimated for these observables and measurements were 
made over 10 5 — 10 6 uncorrelated configurations. 

In the pure Potts case (A = 0), the order parameter 


decays and the susceptibility xs = L 2 ((<S 2 ) — ( S) 2 )/T 
diverges [Fig. [2] close to the transition temperature T c = 
l/log(l + a/3) = 0.9949 [24]. The transition is accom¬ 
panied by a simultaneous proliferation of both types of 
defects, as indicated by an increase in their densities. The 
transition and defect proliferation shift to a higher tem¬ 
perature T ^ 1.26 when vortices are weakly suppressed 
using A = 2 [Fig. [ 2 ]. The location of the transition con¬ 
tinues to shift in this manner, with increasing A, upto a 
threshold value A = A + which we estimate to lie around 
A_|_ ~ 8. 

Above A + , the order parameter clearly exhibits a two 
step decay and the susceptibility shows two distinct peaks 
[Fig-1- For A = 10, the first decay from the ordered 
phase to the intermediate phase occurs around T ^ 1.4 
and is accompanied by the proliferation of domain walls. 
The vortices proliferate near the second decay, which 
marks the transition from the intermediate phase to the 
disordered phase at T « 2.7. Extremely strong suppres¬ 
sion of vortices keeps the first decay unchanged but shifts 
the second decay to T ^ 00 , thus establishing the role of 
vortices in driving the disordering transition. 

Two step decay of magnetization, driven by succes¬ 
sive proliferation of domain walls and vortices, has been 
discussed in the context of Z n vector Potts (clock) spin 
models with n > 5 These models ex¬ 

hibit a phase, intermediate between order and disorder, 
where domain walls proliferate but vortices do not. In 
the intermediate phase, the system fragments into nu¬ 
merous domains in a manner such that the spins fluc¬ 
tuate by arbitrary amounts over large distances and ex- 
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T = 1.2, L = 32 T = 3.0, L = 32 T = 3.0,L = 64 T = 3.0,L = 128 





FIG. 3. (Color online) Top panel shows P(m x ,m y ) in the 
ordered (T = 1.2) and quasi long range ordered (T = 3.0) 
phases for A = 100. Bottom panel shows power-law decay of 
C(r) at different temperatures in the latter phase for L = 256. 
The saturation at large r is due to the finite size of the system. 
77 increases with temperature as shown on the right. 


hibit a U{1) symmetry upon coarse-graining [141. The 
emergent continuous symmetry destroys long range or¬ 
der and gives rise to a quasi long range order which is 
characterized by a power law decay of two-point correla¬ 
tion C(r) = (cos(27r(cro —oy)/3)) oc r -77 , where a 0 and a r 
are spins located at a Eucledian distance r apart on the 
lattice. The exponent 77 changes continuously with tem¬ 
perature throughout the quasi long range ordered phase 
until vortex proliferation disorders the system m ■ 

We have measured the distribution P(m x ,m y ) of the 
Z 3 order parameter m = m x + im yi where m x = 
n a cos( 27 rt/ 3 ) and mn y = n ° sin(27rcr/3) 32 . 34.1. 
The distribution [Fig. [ 3 ] clearly indicates a breaking of 
the three-fold symmetry in the ordered phase and an en¬ 
hancement to U( 1) symmetry in the intermediate phase. 
The U( 1) symmetry survives an increase in system size 
while the magnetization ( m 2 X + m 2 ) 1 / 2 tends to zero in 
accordance with the Mermin-Wagner theorem [35, 155] . 
C(r) exhibits a power-law decay throughout the interme¬ 
diate phase. 77 increases with temperature from 77 ~ 0.35 
and appears to saturate around 77 « 0.75 at high tem¬ 
peratures [Fig. [ 3 ]. This confirms that the intermediate 
phase is indeed a quasi long range ordered phase. 

We now turn to the regime A < 0. The formation of 
vortices is enhanced in this regime and the order-disorder 
transition of the pure Potts case shifts to lower temper¬ 
atures when A is made more negative [Fig. [4]. Addition¬ 
ally, the decay of the order parameter grows sharper and 
the simultaneous rise in the densities of domain walls and 
vortices across the transition becomes more abrupt com¬ 
pared to the pure Potts case. When A is decreased below 
a threshold value A_ ~ —1.3, the model exhibits three 
distinct regimes [Fig. [4]. The densities of domain walls 


and vortices are nearly zero in the ordered phase, show a 
sharp jump to a large value at intermediate temperatures 
and decrease gradually in the disordered phase. The or¬ 
der parameter, on the other hand, shows a sharp decay 
from the ordered phase to the intermediate regime, and 
remains zero thereafter. This might suggest disordered 
behavior in the intermediate regime but an inspection of 
spin configurations in that regime [Fig. [l] reveals that the 
spins are not disordered. Instead, they exhibit a weave 
pattern which corresponds to an ordering of the vortex 
defects: vortices reside on one sublattice and antivortices 
reside on the other. 

In order to characterize the sublattice ordering of the 
vortices, we define a variable e 7 / which is either +1 or -1 
depending on the sublattice of the dual vertex i'. Since 
some of the dual vertices are vacanct (av = 0 ), the sub¬ 
lattice vortex order parameter can be chosen to be of the 
same form as that for a site-diluted Ising antiferromag- 
net 23138]: m svx = L~ 2 JT, gA/ e*/a;*/. The magnitude of 
this order parameter becomes non-zero in the vortex lat¬ 
tice and clearly demarcates the intermediate phase from 
the ordered and disordered phases [Fig. [ 4 ]. Since this 
order parameter, based on topological defects, is able to 
distinguish the vortex lattice phase from the disordered 
phase, while the symmetry based order parameter S is 
unable to do so, the vortex lattice phase provides a simple 
example of classical topological order. Formation of vor¬ 
tex lattices in superfluids and superconductors has been 
extensively studied over the past few decades HESHIHin 
superconducting thin films, proliferation of dislocations 
and disclinations drive a two step structural melting of 
the vortex lattice mmm- In the present model, the 
melting of the vortex lattice to the disordered phase oc¬ 
curs via a single step process as indicated by the decay of 
(\m svx \) [Fig. 0 - The sublimation of the vortex lattice to 
the ordered phase, which occurs at a lower temperature, 
is marked by a sharp decay of (\m svx \). With a further 
decrease of A, the melting shifts to higher temperatures 
[Fig- E]. The sublimation, on the other hand, shifts to 
lower temperatures and hits the zero temperature limit 
at A « —1.5. Below this A, the ordered phase is absent 
and the model exhibits a single melting transition from 
the vortex lattice to the disordered phase. 

The nature of the new phases uncovered above is quite 
clear from the data obtained for small systems. A pre¬ 
cise estimate of the temperature and nature of the tran¬ 
sitions, on the other hand, requires detailed analysis of 
data from large systems and will be presented separately. 
Here, we make a few comments regarding the possible 
nature of the transitions. The order-disorder transition 
for A = 0 is of the second order type [24] [45]. This 
order-disorder transition continues to occur throughout 
the range A_ < A < A + and is accompanied by a cou¬ 
pled proliferation of domain walls and vortices. However, 
the decay of the order parameter appears to grow weaker 
with increasing A [Fig. [ 2 ] and sharper with decreasing 
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FIG. 4. (Color online) Weak enhancement of vortices (A = — 1) results in a sharp decay of S accompanied by simultaneous 
increase of pdw and p vx . Stronger enhancement (A = —1.4) opens up an intermediate vortex lattice phase characterized by a 
non-zero value of (\m svx \}- With further enhancement, the ordered phase vanishes while the vortex lattice melts at a higher 
temperature. Data corresponds to system sizes mentioned in Fig. [2] 


A [Fig.g. This observation leads us to conjecture that 
the critical exponents of the transition might vary with 
A. The four state (Ashkin-Teller) ferromagnet is known 
to exhibit continuously varying criticality along a line 
of order-disorder transitions due to an interplay between 
vortices and domain walls [46]. The effect of the coupled 
proliferation on the nature of the transition in the present 
model, therefore, seems to be an interesting problem. In 
this context, we note that claims of continuously varying 
criticality in the closely related three state chiral Potts 
model continues to be a controversial issue [474149] . 

For A > A + , the order-disorder transition splits into 
two and the intermediate phase exhibits quasi long range 
order. In Z n ferromagnets with n > 5, the two transitions 
bordering the quasi long range ordered phase are known 
to be of the Berezinskii-Kosterlitz-Thouless type mm- 
We can expect the two transitions in the present model to 
be of that type as well. In Z n models, the decay exponent 
r] changes from rj = 4/ n 2 at the low temperature transi¬ 
tion to r] = 1/4 at the high temperature transition [34] . 
The values of p E (0.35,0.75) obtained in the quasi long 
range ordered phase of the present model [Fig. [3] fall be¬ 
yond that range. This deviation from the standard Z n 
model scenario indicates that the bounds for p can change 
with vortex suppression. 

Our estimate of A + « 8 is an approximate one. For A = 
6, the intermediate phase is narrow in small systems and 
appears to shrink with increasing L [Fig. [2]. On the other 
hand, for fixed L, the extent of the phase increases with 
with increasing A. This competing effect of L and A offers 
two possibilites: (a) there exists a threshold A + , above 
which the intermediate region has a non-zero extent in 
the thermodynamic limit, or (b) the intermediate region 
shrinks to a point in the thermodynamic limit for all 


finite A and the model exhibits an extended quasi long 
range ordered phase only in the A —>• oc limit. A precise 
estimate of A + , therefore, remains an open problem and is 
reminiscent of the problem regarding the location of the 
Lifshitz point in the three state chiral Potts model [Sol- 
154]. 

In the range A_ < A < 0, the order-disorder transition 
grows sharper with decreasing A, hinting at the possibil¬ 
ity that the transition will become discontinuous before 
A goes below A_. Such a scenario is quite plausible be¬ 
cause an abrupt proliferation of vortices, similar to the 
behavior shown in Fig. [4] is known to induce discontin¬ 
uous behavior H3E5]. For A < A_, the transition from 
the ordered phase to the vortex lattice phase [Fig. [I] is 
clearly discontinuous. The gradual decay of (\m svx \) be¬ 
tween the vortex lattice phase and the disordered phase, 
on the other hand, suggests that the melting transition 
might be second order in nature. 

The universality class of the three state Potts tran¬ 
sition has a ubiquitous presence in the physics of sta¬ 
tistical [24] [564164] . quantum [65 69 and gauge sys¬ 
tems |31 E0HZ1- We have shown that the transition 
is driven by a coupled proliferation of domain walls and 
vortices. By manipulating the formation of the defects, 
we have uncovered two new phases in the model. Apart 
from the exciting possibility that these phases might be 
realizable in some of the systems, our work provides a 
step towards better understanding the role of topological 
defects and the presence of topolgical order in classical 
spin models. 
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